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ABSTRACT 

Aims. We aim at characterizing a sample of nine new hard X-ray selected Cataclysmic Variable (CVs), to unambiguously identify 
them as magnetic systems of the Intermediate Polar (IP) type. 

Methods. We performed detailed timing and spectral analysis by using X-ray, and simultaneous UV and optical data collected by 
XMM-Newton, complemented with hard X-ray data provided by INTEGRAL and Swift. The pulse arrival time were used to estimate 
the orbital periods. The broad band X-ray spectra were fitted using composite models consisting of different absorbing columns and 
emission components. 

Results. Strong X-ray pulses at the White Dwarf (WD) spin period are detected and found to decrease with energy. Most sources are 
spin-dominated systems in the X-rays, though four are beat dominated at optical wavelengths. We estimated the orbital period in all 
system (except for IGR J16500-3307), providing the first estimate for IGR J08390-4833, IGRJ18308-1232, and IGRJ18173-2509. 
All X-ray spectra are multi-temperature. V2069 Cyg and RX J0636+3535 posses a soft X-ray optically thick component at kT~80 eV. 
An intense 1% Fe line at 6.4 keV is detected in all sources. An absorption edge at 0.76 keV from OVII is detected in IGR J08390-4833. 
The WD masses and lower limits to the accretion rates are also estimated. 

Conclusions. We found all sources to be IPs. IGR J08390-4833, V2069 Cyg, and IGR J16500-3307 are pure disc accretors, while 
IGR J 1 8308- 1232, IGR J1509-6649, IGR J 1 7 195-4100, and RX J0636+3535 display a disc-overflow accretion mode. All sources 
show a temperature gradient in the post-shock regions and a highly absorbed emission from material located in the pre-shock flow 
which is also responsible for the X-ray pulsations. Reflection at the WD surface is likely the origin of the fluorescent iron line. There 
is an increasing evidence for the presence of a warm absorber in IPs, a feature that needs future exploration. The addition of two 
systems to the subgroup of soft X-ray IPs confirms a relatively large (~ 30%) incidence. 

Key words. Stars: binaries: close - Stars: individual: IGR J08390-4833, IGR J18308-1232, IGR J16500-3307, IGR J18173-2509, 
IGR J17195-4100, V2069 Cyg (also known as RX J2 1237+42 1 8), RX J0636+3535 (also known as V647 Aur), IGR J15094-6649, 
and XSS J0056+4548 (also known as V515 And) - X-rays: binaries - Accretion, accretion discs - starsmovae, cataclysmic variables 


1. Introduction 

Magnetic Cataclysmic Variables (CVs) constitute a subgroup of 
the CV class, harboring accreting white dwarfs (WD) with mag- 
netic field strengths B S 10 5 G. These systems are further sub- 
divided in two groups, depending on the WD magnetic field in- 
tensity and degree of synchronism (P ro t=a>/Porb=n)- Those called 


* Based on observations obtained with XMM-Newton and 
INTEGRAL, ESA science missions with instruments and contri- 
butions directly funded by ESA Member States; and Swift, a NASA 
science mission with Italian participation. 


Polars are synchronous mCVs. They show signatures of strong 
magnetic fields (B ~ 10 - 230 MG) through the presence of 
a conspicuous polarization at optical and near-IR wavelengths. 
The so-called Intermediate Polars (IPs) possess instead asyn- 
chronously rotating WDs (P w << Pq) and except for a few cases 
(Piirola et al. 1993; Potter et al. 1997; Katajainen et al. 2007; 
Butters et al. 2009; Potter et al. 2012), they do not show polar- 
ization. These systems are believed to possess weakly magne- 
tized WDs (B <10 MG). The Polars populates the orbital period 
distribution at short periods mainly below the 2-3 h orbital pe- 
riod ’’gap” (Warner 1995; Wheatley 1995). On the other hand. 
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IPs are generally found above the gap, with only a handful of 
systems below it. 

The material lost from the Roche-lobe overflowing late- 
type secondary star (typically a Main Sequence or a sub-Giant) 
may flow toward the WD in different modes, depending on the 
magnetic field intensity: directly (stream), through a truncated 
disc or ring. The high field Polars are direct accretors, while 
IPs may accrete via a stream or a disc/ring depending on their 
magnetic moment and spin-to-orbital period ratio (Norton et al. 
2004, 2008). A combination of these, called disc-overflow, may 
also occur (Hellier 1995; Norton et al. 1997) and is linked to the 
mass accretion rate. Close to the WD, the magnetically con- 
fined flow of matter is accreted through a column in stream- 
fed systems and through an arc-shaped curtain (Rosen et al. 
1988) in disc/ring-fed systems. Approaching the WD surface, 
the velocities are supersonic and a stand-off shock is formed. 
In the post-shock region (PSR) the gas at a temperature of 
kT~ 10-60 keV slows and cools via thermal bremsstrahlung and 
cyclotron radiations (Aizu 1973; Wu et al. 1994; Cropper etal. 
1999). The efficiency of the two cooling mechanisms depends 
on the WD magnetic field strength (Woelk & Beuermann 1996; 
Fischer & Beuermann 2001). The post-shock emission is par- 
tially intercepted and thennalized by the WD surface, giving rise 
to an optically thick emission that emerges in the soft X-rays 
and EUY/UV regimes. An intense soft (kT ~ 20-50 eV) com- 
ponent is a characteristic signature of the X-ray spectra of the 
Polars, mostly balancing the cyclotron flux (Beuermann 2004). 
On the other hand IPs have typically stronger hard X-ray fluxes 
than the Polars and only recently a soft highly absorbed opti- 
cally thick component has been detected (de Martino et al. 2004; 
Evans & Hellier 2007; Anzolin et al. 2008, 2009). 

Though mCVs represent a relatively small (~ 20%) frac- 
tion of CVs, this number is rapidly increasing thanks to the re- 
cent hard X-ray surveys conducted by INTEGRAL and Swift 
above 20keV. So far 64 CVs are identified in the latest cat- 
alogue releases (Bird etal. 2010; Cusumano et al. 2010), 43 
of them are mCVs. The IPs represent ~ 80% of this hard 
X-ray sample (including both already known (Barlow et al. 
2006) and recently discovered members (Bonnet-Bidaud et al. 
2007; de Martino et al. 2008; Anzolin et al. 2009; Butters et al. 
2008; Bonnet-Bidaud et al. 2009; Pretorius 2009; Scaringi et al. 
2011), though one misidentification was found (de Martino et al. 
2010). Hence, while in the pre-INTEGRAL and Swift times IPs 
amounted to only ~20 members, this number has more than dou- 
bled nowadays. 

Hard X-ray mCVs have the potential to be important con- 
tributors to the X-ray source population at low luminosities 
(~ 10 30 - 10 33 erg s -1 ). They were proposed to be the major con- 
stituent of galactic ridge (Sazonov et al. 2006; Revnivtsev et al. 
2008, 2009) and galactic bulge (Revnivtsev et al. 2011; 
Hong et al. 2012) X-ray emission. Therefore, they are also be- 
lieved to have an important role in the X-ray luminosity function 
of other galaxies. 

In the framework of a program conducted with 
XMM-Newton, with the main goal to identify new mag- 
netic type CVs, we present here the analysis of simultane- 
ous X-ray, UV and optical data of a sample of nine CVs 
that were detected as hard X-ray sources by INTEGRAL, 
RXTE or SWIFT: IGR J08390-4833, IGR J18308-1232, 
IGR J1 8 1 73-2509, IGR J17195-4100, V2069 Cyg (also 
known as RX J2 1237+42 18), V647 Aur (also known as 
RX J0636+3535), IGR J15094-6649, V515 And (also known as 
XSS J0056+4548), and IGR J16500-3307 (hereafter IGR J0839, 
IGR J1830, IGR J1817, IGR J1719, V2069 Cyg, RX J0636, 


IGR J1509, XSS J0056, and IGR J1650 respectively. We com- 
plement the X-ray analysis with a high-energy coverage from 
INTEGRAL/IBIS and Swift/BAT publicly available data. Five 
of these sources (IGR J0839, IGR J1 8 17, XSS J0056, IGR J 1509 
and IGR J1 71 9) have also Chandra, RXTE, or Swift/XRT 
coverage. All, except IGRJ0839, IGRJ1830 and IGRJ1817, 
were also observed in optical photometry and spectroscopy 
from ground-based telescopes, from which orbital periods 
were detennined and optical pulses detected, thus allowing a 
comparison with our results. 

In Sect. 3.1 we present the X-ray, UV and optical timing 
analysis. In Sect. 3.2 we present the analysis of the X-ray broad- 
band spectra of each source. For the three brightest sources 
we also present the high resolution spectra provided by the 
XMM-Newton RGS instrument. In Sect. 4 we discuss the ori- 
gin of rotational pulses, the accretion mode and the emission 
properties of our sample. 

2. Observations and data analysis 

2.1. XMM-Newton observations 

The XMM-Newton Observatory includes three ~ 1500 cm 2 
X-ray telescopes with an EPIC instrument in each focus, 
a Reflecting Grating Spectrometer, RGS (den Herder et al. 
2001) and an Optical Monitor, OM (Mason et al. 2001). 
Two of the EPIC imaging spectrometers use MOS CCDs 
(Turner etal. 2001) and one uses a PN CCD (Striider et al. 
2001). XMM-Newton collected data of all nine CVs. The main 
observation parameters, together with that of all other observa- 
tory instruments are reported in table 1 . 

Data were processed with SAS version 10.0.0, using the 
updated calibration files (CCF) available in January 2011. All 
observations were performed with the EPIC-PN (PN hereafter) 
camera set in prime full window imaging mode (time resolu- 
tion^. 0734 s) with a thin filter applied. Standard data screening 
criteria were applied in the extraction of scientific products. We 
accumulated a one-dimensional image and fitted the ID photon 
distribution with a Gaussian. Then, we extracted the source pho- 
tons from a circular region of radius 40" centered at the Gaussian 
centroid. The background was obtained (within the same CCD 
where the source lies) from a circular region of the same size. 
Single and double pixel events with a zero quality flag were se- 
lected for the PN data. For the spectral analysis we cleaned all 
observations from solar flares by collecting CCD light curves 
above 10 keV and applying an intensity threshold. Spectra were 
then produced only for those parts of observations under the 
threshold limit. On the other hand, for the timing analysis we 
used the epiclccor task, producing a background subtracted 
light curve in the range 0.3-15 keV (with a bin time of 1 1 s), con- 
sequently removing contamination from solar flares. The source 
event arrival time of each observation, in the 0.3-15 keV energy 
range, were converted into barycentric dynamical times (BDT) 
by means of the SAS tool barycen. 

The PN spectra were rebinned before fitting, to have at least 
30 counts per bin. We report the analyses obtained with the PN 
data only (consistence with the results of EPIC -MOS cameras 
was always verified). Phase-resolved spectra were also extracted 
at the pulse maximum and minimum. All spectra were analyzed 
using the version of XS PEC (12.5. In). 

We also inspected the RGS spectra. We used the 1st order 
spectra and responses produced by the pipeline, which runs an 
automatic version of the RGS meta-task, rgsproc. While, in 
principle, the quality of data products can be improved some- 
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what by manually repeating the reduction, we found no obvious 
problems with the pipeline products. Also considering the lim- 
ited statistical quality of the RGS data for our sample, we believe 
our choice to use the standard pipeline products is well justified. 
The RGS spectra ofIGR J1719, IGR J1509, andXSS J0056 were 
of high enough statistical quality to merit further analysis; these 
are discussed below in sect. 3.2.3. The data for the other targets 
were of too low S/N to add useful information. We grouped the 
RGS spectra of the three sources such that each new bin had at 
least 25 counts. 

Simultaneous coverage to the X-rays was ensured by the 
OM operated in fast window mode using sequentially the B 
(3815-4910A) or the U (3000-3800 A) filter and UVM2 (2000- 
2800 A) or UVW1 (2450-3200 A) filter. Consequently, each tar- 
get was observed for half of the time of the EPIC exposure in 
the B or U filter and the other half in a UV filter. Due to the 
OM fast window mode acquisition, the OM photometric data 
in each filter consist of sequential (four or five) segments of 
length ~ 1300-2800 s each (according to the total EPIC expo- 
sures). Average net count rates and instrumental magnitudes are 
reported for those sources with usable data in table 1, where we 
also report whether no variability is detected. Background sub- 
tracted OM light curves were obtained using the standard SAS 
processing pipeline with a binning time of 1 0 s. Barycentric cor- 
rections were also applied to the OM light curves. 


3.1.1. X-ray variability 

We computed the power spectra in the 0.3-15 keV range using 
the PN data. In all of them, strong peaks are found and in many 
cases also harmonics (see figure 1, 2, and 3). The periods of the 
main peaks were detennined by means of a phase-fitting tech- 
nique, see Dall’Osso et al. (2003) for details on the technique. 
The derived values are reported in table 2. All uncertainties are 
hereafter at lcr confidence level (c.l.) if not otherwise specified. 
For IGR J1 71 9 we did not detect significant periodic variability 
in the last third of the observation and therefore, the results refer 
to the first 2/3 only. This source also shows a variability on time 
scale of a few thousands of seconds occurring at about 1/3 of 
the total observation time. This variability is likely responsible 
for the low frequency peaks seen in the power spectrum. We are 
however unable to characterize it with the present data. 


2.2. The INTEGRAL observations 

The INTEGRAL IBIS/ISGRI instrument (Ubertini et al. 2003; 
Lebrun et al. 2003) observed all sources except RX J0636. Hard 
(20-100 keV) X-ray data were extracted from all pointings 
within 12° from the source positions. 

To study the persistent X-ray emission, the time-averaged 
ISGRI spectra were obtained from mosaic images in five energy 
bands, logarithmically spaced between 20 and 100 keV. Data 
were reduced with the standard OSA software version 7.0 and 
then analyzed with the algorithms described by Goldwurm et al. 
(2003). 

2.3. The Swift observations 

The Swift Burst Alert Telescope, BAT (Barthelmy 2000), is a 
wide-field (~1 steradian) coded aperture mask instrument sensi- 
tive in the 14-195 keV range. Thanks to the large field of view, 
BAT has built up a sensitive all-sky map of the hard X-ray sky. 
The BAT team provides the average 8-channel spectrum of each 
source detected over the first 58 month of the mission. 

RX J0636 only recently appeared in the Swift/BAT 
58 month survey catalogue. To obtain an high- 

energy coverage of this source we downloaded 
the BAT spectrum along with the response matrix 
(http://swift.gsfc.nasa.gov/docs/swift/results/bs58mon/). 


3. Results 

3.1. Timing analysis 

All sources show periodic-like variability in their EPIC light 
curves and six of them also in their OM light curves. We then 
analyzed these light curves to identify periodic signals and to 
characterize them at different energies. 
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Table 1 . Summary of main observations parameters for all instruments. 


Source 

Telescope 

OBSID 

Instrument* 

Date 

yyyy-mmm-dd 

UT s tart 

hh:mm 

T 

L expo 

(ks)** 

Net count rate (mag)*** 
cts/s 

IGR J08390-4833 

XMM-Newton 

0651540101 

EPIC-pn 

2010 Dec 20 

12:34 

33.9 (27.3) 

0.776 ± 0.006 




OM-B 


12:20 

15.0 

3.6 + 0.9(17.8) 




OM-UVM2 


16:57 

15.2 

0.25 + 0.05 (17.3) 


INTEGRAL 


IBIS/ISGRI 



2000 

0.26 + 0.04 

IGR J 1 8308-1232 

XMM-Newton 

0601270501 

EPIC-pn 

2010 Mar 11 

15:26 

27.6(15.7) 

0.793 + 0.007 


INTEGRAL 


IBIS/ISGRI 



2200 

0.36 ± 0.05 

IGR J16500-3307 

XMM-Newton 

0601270401 

EPIC-pn 

2010 Feb 2 

11:04 

30.0(18.0) 

1.33 + 0.01 




OM-B' 1 


10:50 

6.2 

10 + 2(16.7) 




OM-UVW1" 


13:00 

8.0 

2.2 + 0.5 (16.3) 


INTEGRAL 


IBIS/ISGRI 



2200 

0.31 +0.02 

IGR J1 8 1 73-2509 

XMM-Newton 

0601270301 

EPIC-pn 

2009 Sep 07 

02:03 

35.0 (30.0) 

0.941 + 0.006 




OM-B 


01:49 

15.7 

3.0+ 1.0(18.1) 


INTEGRAL 


IBIS/ISGRI 



4500 

0.88 + 0.04 

IGR J1 7 1 95-4 1 00 

XMM-Newton 

0601270201 

EPIC-pn 

2009 Sep 03 

06:58 

24.0" (27.0) 

6.70 + 0.02 




OM-B 


06:44 

11.4 

36 + 5(15.4) 




OM-UVM2 4 


06:44 

11.4 

2.3 + 0.6(14.9) 




RGS 1&2 


06:35 

33.9 



INTEGRAL 


IBIS/ISGRI 



2600 

0.93 + 0.05 

V2069 Cyg 

XMM-Newton 

0601270101 

EPIC-pn 

2009 Apr 30 

11:09 

26.4(12.4) 

1.05+0.01 




OM-B 


10:55 

10.6 

10+1.5 (16.8) 




OM-UVM2" 


14:18 

10.6 

0.17 + 0.04(17.6) 


INTEGRAL 


IBIS/ISGRI 



1900 

0.14 + 0.05 

RX J0636+3535 

XMM-Newton 

0551430601 

EPIC-pn 

2009 Mar 18 

17:01 

29.2(12.1) 

1.37 + 0.01 




OM-B" 


16:47 

6.7 

13 + 2(16.4) 




OM-UVM2" 


22:41 

6.7 

1.5 + 0.2(15.4) 


Swift 


BAT 



~ 7500 

0.00028 ± 0.00003 

IGR J1 5094-6649 

XMM-Newton 

0551430301 

EPIC-pn 

2009 Feb 02 

13:39 

30.0 

2.56+0.01 




OM-U 


13:25 

14.2 

36 + 4(14.3) 




OM-UVM2" 


17:49 

12.4 

2.8 + 0.8(14.7) 




RGS 1&2 


13:16 

31.9 



INTEGRAL 


IBIS/ISGRI 



1400 

0.61 ±0.06 

XSS J0056+4548 

XMM-Newton 

0501230301 

EPIC-pn 

2007 Dec 31 

04:56 

15.0 (7.4) 

3.06 ± 0.02 




OM-B 


04:42 

6.7 

69 + 5(14.6) 




OM-UVM2 


07:01 

6.7 

10+ 1 (13.2) 




RGS 1&2 


04:33 

16.9 



INTEGRAL 


IBIS/ISGRI 



240 

0.35 ± 0.06 


* Optical Monitor band reported when data usable 

** Net exposure times. In parenthesis is reported the solar flare removed exposure 

*** Instrumental magnitude reported in parenthesis. 

" No variability detected. 

4 Time window filter applied: only first 2/3 of observation are used. 
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Table 2. Timing properties of the source sample. From left to right: (X-ray main (spin) period); Pf ide (X-ray sidebands); P^ s,de (orbital period derived from X-ray sidebands); 

(orbital period derived from phase-fitting); A^/A^ e (spin to sideband amplitude ratio); optical spin period from OM; P’^ M (optical beat from OM); p" pl ’' per 

(optical spectroscopic period with its reference in parenthesis); P'^ (adopted orbital period in this work. 
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XSS J0056+4548 



Time since BJD=2454465.5 [days] 

Fig. 4. PF vs time for XSS J0056. The sinusoid period with a 30- 
uncertainty is also reported. 


For each source, the light curve in the 0.3-15 keV range was 
folded at the main (P w ) period reported in table 2, and then fitted 
with a series of sinusoids (one for each harmonic). The statistical 
significance of the inclusion of higher harmonics with respect 
to the fundamental one was evaluated by an F-test (see table 
3). The first harmonic was always statistically significant (here- 
after we consider significant a probability greater then 3 cr c.l.), 
with the exception of IGRJ1830 and IGRJ1719 in which only 
the fundamental is significant. The second harmonic (3 at) was 
found to be significant in IGR J0839, V2069 Cyg and IGR J1 509. 
V2069 Cyg is also the only source showing 4a>. 

The power spectra of most sources either show an excess of 
power around the fundamental frequency or additional weaker 
peaks close-by. In some cases also the harmonics show similar 
excess or peaks. Therefore, to estimate their frequency we de- 
trended the light curves from the main pulse by subtracting the 
composite sinusoidal function corresponding to the fundamen- 
tal frequency plus their harmonics when detected. The residual 
power spectra of the detrended light curves makes clearly visi- 
ble additional peaks. Consequently, we fitted the detrended light 
curves with one new sinusoid (or more sinusoids when harmon- 
ics are present). We reported the estimated period, when they re- 
sulted to be statistically significant, in table 2. These additional 
periodicities are identified as the orbital sidebands of the main 
pulse (at). 

IGR J1 7 19, RXJ0636 show the negative sideband a> - £1 
(beat). In IGR J1 509 we detected the harmonic of the beat 
3 a> - £1, but not the fundamental, and the 2a> - £1 sideband. In 
RXJ0636 we also detected at - 2£1. For this source the pres- 
ence of both negative sidebands is significant at 4<x c.l.. For 
IGRJ1830 the at - 2£i is detected at 4.1<x c.l., but the at - £1 
is only at 1.5<x c.l.. The positive at + £2 sideband is detected in 
XSS J0056 (7.4o-c.l.) and2(w+ST) in IGR J1817 (16o-c.l.) since 
the main X-ray period is the first harmonic (see sect. 3.1.3). 

We here note that Sazonov et al. (2008) using a short 
Chandra observation of IGRJ0839 detected a periodicity of 
1450 ± 40 s consistent with our more constrained value. 

Butters et al. (2008) failed to detect a periodicity in a RXTE 
observation of IGR J1 7 19. On the other hand, they detected 


two periodic signals in XSSJ0056 again using RXTE. One at 
465.68±0.07 s and a weaker one at 489.0±0.7 s. The former is 
inconsistent with our determination (being out by more than 5 
s), while the latter is not detected in our data. 

For IGR J1 8 1 7 Nichelli et al. (2009) using a Swift/XRT pointing 
report a period of 830.7±0.07 s consistent with our results and a 
period of ~1660 s from a previous short Chandra exposure (see 
sect. 3.1.3 for its interpretation). 

Butters etal. (2009) by using RXTE data of IGR J 1509 
found a signal at 809.7±0.6 s, consistent with our results. 

The phase-fitting procedure also allowed us to get an esti- 
mate of the orbital period. The light curve of each source was 
subdivided in n intervals (10, 15 or 20), where n depends on 
the source brightness and on the intensity of the signal. Then 
we folded each light curve at the period reported in table 
2 and then we fitted it with a Fourier sine series truncated at 
the highest significant harmonic. We consequently studied the 
variations with time of the phase (and amplitude) of the main 
pulsation, corresponding to the fundamental frequency. In most 
cases the phase variations, typically within A0± 0.2, appear to 
be randomly distributed around the average value. Only in three 
sources, IGR J1 8 17, RX J0636 and XSS J0056 the lag variations 
show a clear trend, allowing a sinusoidal fit. These are shown in 
figure 1, 2, and 3. The statistical significance of the inclusion of 
a sinusoid respect to a constant function is evaluated with an F- 
test, and was found to be 3.41, 3.00 and 4.21<x respectively. The 
resulting periods are reported in column 5 of table 2. 

Most sources show also variations with time of the pulse am- 
plitude (or PF, see sect. 3.1.2 for its definition), but these ap- 
pear to be randomly distributed around an average value, with- 
out showing a clear trend. An exception is XSSJ0056, where 
also the amplitude changes following a sinusoidal modulation. 
A sinusoidal fit is statistically significant at 3.52<x (see figure 4), 
giving a period of 2. 6 ±0.4 h (3<xc.l.). Hence, while in RX J0636 
and IGR J1817 the phase shifts could be due to time arrival de- 
lays, in XSS J0056 the variations seen in both amplitudes and 
phases are likely intrisic. This is also suggested by the presence 
of a peak around 9.3 x 10 -5 Hz in the power spectrum and a si- 
nusoidal fit gives a period of 2.52±0.08 h, in agreement with the 
spectroscopic period found by Bonnet-Bidaud et al. (2009). 

Interpreting the main pulsation as the WD spin period, the 
sidebands can be used to estimate the binary orbital periods, re- 
ported in column 4 of table 2. For IGRJ1719 and IGRJ1509 
the periods are in agreement with those derived from spec- 
troscopy, reported for comparison in column 9. The values found 
for IGRJ1817 and XSS J0056 are consistent within errors with 
those derived from the phase-fitting method, but for XSS J0056 
they are not compatible with the results from the PF vs time anal- 
ysis and from the sinusoidal fit. For RXJ0636 we found con- 
sistency between the period found from phase-fitting and that 
derived using the beat but not at - 2£1. However, all derived 
values are largely off from the tentative period of 3.4 h found 
from optical spectroscopy (Gansicke et al. 2005). Furthermore, 
V2069 Cyg also shows in its power spectrum substantial power 
at low frequencies, but we are unable to constrain a period from 
the light curve. However, a sinusoidal fit with period fixed at 7.5 
h, derived from spectroscopy by Thorstensen & Taylor (2001), 
provides a possible match with data and gives a 35% amplitude 
modulation. This variability is more pronounced at soft energies 
(< 3 keV), where we also detect a dip, lasting ~1500 s, around 
the minimum of this long-term X-ray modulation (figure 5). 

For the systems under study we hereafter adopt, when avail- 
able, the more accurate orbital spectroscopic periods, except for 
RXJ0636 for which the orbital period is likely larger, Pq <;5 
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IGR J0839-4833 


IGR J0839— 4833 



Frequency (Hz) 
IGR J1830-I232 





J 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 L 
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Time since BJD=2455550.5 [days] 

IGR J1 830-1232 




0.7 0.8 0.9 

Time since BJD=2455266.5 [days] 


IGR J1650— 3307 


IGR J1 650-3307 




Fig. 1 . Left Panel : Power spectrum. Fundamental and harmonics are reported in red, while dotted blue lines represent sidebands. 
Right panel : Evolution of the phase of the main signal, a>, vs time. 
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Fig. 2. As in figure 1 . P or b with 1 cr uncertainty is also reported for IGR J 1 8 1 7 . 
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Fig. 3. As in figure 1 . For IGR J1509, the two dotted blue lines refer to 2a> - Q and 3 (a> - Cl) respectively. P or b is also reported for 
RX J0636 (with 3cr uncertainty) and for XSS J0056. 
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Table 3. Pulsed fraction vs energy. Highest detected hannonic and its statistical significance is also reported. 


Source 



Pulsed Fraction 


not 

sign. 


0.3-1 keV 

1-3 keV 

3-5 keV 

5-15 keV 

0.3-15 keV 




% 

% 

% 

% 

% 


CT 

IGR J0839 

99 ±2 

74 ± 1 

35 + 2 

18 + 2 

57+ 1 

3 

3.97 

IGR J 1830 

30 ±4 

33+2 

26 ±2 

21+2 

28+ 1 

1 

- 

IGR J 1650 

18 ± 2 

12 ± 1 

5 ± 2 

6 + 2 

9+1 

2 

4.29 

IGR J1817" 

58 + 3 

46 ±2 

34 ±2 

27+1 

30+ 1 

2 

5.94 

IGRJ1719 6 

14 ± 1 

13 ± 1 

6+ 1 

< 3 

10+ 1 

1 

- 

V2069 Cyg 

14 + 3 

25 ±2 

16+ 1 

9+ 1 

18+ 1 

4 

3.5 

RX J0636 

18 + 2 

21 ±2 

< 

17 

15 ± 1 

2 

3.02 

IGR J 1509 

57 ±1 

51 ± 1 

27+1 

9+ 1 

42+ 1 

3 

3.52 

XSS J0056 

42 + 2 

31 ± 1 

16 ± 2 

9 + 1 

24+ 1 

2 

4.54 


“ Pulsed fraction refers to Id. 
b Derived from the first 2/3 of observation. 
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Fig. 5. 0.3-1 keV lightcurve of V2069 Cyg binned at 198 s. 
Reference time is the PN start of exposure, see table 1. A dip 
is present at about 1.5 x 10 4 s. 


h. For those sources without a spectroscopic determination we 
adopt the orbital periods as derived from the X-ray sidebands 
(for IGR J0839 from optical sideband). Exception is IGR J1817, 
for which we adopt the range 6. 3-8. 7 h, resulting from sidebands 
and phase-fitting method. The adopted values are reported in the 
last column of table 2. 

3.1.2. Pulsed fraction 

We also studied the energy dependence of pulses. We selected 
five energy bands, producing four background subtracted light 
curves binned at 66 s (0.3-1, 1-3, 3-5, 5-15 keV), plus the 0.3- 
15 keV light curve with a bin time of 1 1 s, and folded them 
at the value of P* reported in table 2. For all sources the X- 
ray pulse profiles have structured pulse shapes. Moreover, for all 
sources the pulse profiles resulted to be roughly phase aligned at 
all energies, within statistical uncertainties (see figure 6 and 7). 


The folded light curves were then fitted with a Fourier sine 
series truncated at the highest significant harmonic. This allowed 
to study, for each energy interval, the pulsed fraction (PF), here 
defined as: PF = (A max - A min )/(A max + A min ), where A max and 
A m i n are respectively the maximum and minimum value of the 
sinusoid used to account for the fundamental harmonic. For all 
sources the PF was found to clearly decrease with energy (see 
table 3). In order to further inspect for possible spectral changes 
with the spin phase, we also construct the Hardness Ratios (HRs) 
using the three energy intervals 0.3-1 keV 1-3 and 3-5 keV. 
The HR is here defined as: HR = n l te v/ n I ce vb’ where nj^y is the 
number of photons detected by the PN camera in the inspected 
energy range (called a and b), and i is the phase interval. From 
HR = n‘ 15 /n' 13 we found for all sources a hardening at the 
pulse minimum, except for IGR J1 830 and V2069 Cyg, where 
HR was found to be constant. For V2069 Cyg, due to its peculiar 
spectrum in the soft energy range (E< 3 keV; see sect. 3.2.1), we 
also inspected the HR = n 1 1 _ 3 /nj )3 _ 1 , that indicate a softening at 
spin minimum, suggesting the contribution of a soft component. 

3.1.3. The Optical and UV variability 

Optical and UV light curves were inspected for the presence of 
periodic signals by producing power spectra in each band, except 
for RX J0636, IGR J1830 and IGR J1650 that were too faint. The 
UV light curves are generally noisier and for most of them we 
are unable to find significant periodicities. 

The light curves were also fitted with a sinusoidal function. 
The periods were always found to be consistent with one of those 
detected in the X-rays (see table 2). Worth noticing are the fol- 
lowing cases: 

IGR J0839: The UV light is modulated (PF=35%) at a period 
(1480 ± 13 s) consistent with that found in the X-rays, but the 
B band light curve shows a period of 1560 ± 7 s (PF=22%). 
The latter is likely the orbital sideband (beat) a>-Q, suggesting an 
orbital period of Pn=8 ± 1 h. Hence this source is spin-dominated 
in the X-rays and UV, but beat-dominated in the optical. 

IGR J1 8 17: From the light curve in the B filter we derived a 
period ofl690±18s, twice the detected X-ray period. Moreover, 
no signal was found at the X-ray period. We conclude that the 
X-ray period corresponds to the first harmonic, while the optical 
period is the spin. 
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Fig. 6. X-ray pulse profiles at different energy intervals for different sources. Energy increases from top to bottom. Optical and UV 
band pulse profiles are also reported when available. 
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Fig. 7. As in figure 6. 


XSSJ0056 + 4548 



Spin Phase 


Fig. 8. OM light curve of XSS J0056 in the B band for the five 
time intervals of exposure. Time increases from bottom to top. 


IGRJ1719: The B band light curve, simultaneous with the 
X-rays, shows a period of 1 102 ± 4 s. This value falls in be- 
tween the two periods detected in the X-rays (1062 ± 2 s and 
1 163 ± 4 s), suggesting that both periodicities are also present at 
optical wavelengths. This is also confirmed by ground-based op- 
tical photometry (Pretorius 2009) that showed a weaker ~ 1055 
s and a stronger ~ 1 139 s peak in the power spectra. Therefore, 
also this system is spin-dominated in the X-ray, while beat- 
dominated at optical wavelengths (see also sect. 4.2). 

V2069 Cyg: we derived a pulse period of 751 ± 5 s, consis- 
tent, within errors, with that found in the X-rays. Interestingly, 
folding the light curve at the X-ray period, the optical pulse is 
anti-phased with the X-ray one. Moreover, it shows a broad max- 
imum centered on the X-ray minimum, with a dip where an X- 
ray secondary maximum occurs (see Fig 6). 

XSS J0056: it has peculiar UV and optical variabilities. The 
power spectra in these bands did not show significant peaks. 
However, inspection of the individual five segments (~1.34ks 
each) in the two bands revealed the presence of non-stationary 
pulses. They are found to be phase aligned with the X-ray pulse 
during the first 2/5 of the OM B band exposure, for which a sinu- 
soidal fit gives a period of 474±4s (consistent with what found 
in X-rays), but antiphased in the third and fifth segment (see fig- 
ure 7 and 8). The UVM2 filter data revealed a weak variability 
(~ 9%) during the third and possibly the fourth segment, but we 
are unable to determine the period. 

3.2. Spectral analysis 
3.2.1. 0.3-100 keV spectrum 

For all sources, a single-temperature model is inadequate to de- 
scribe the broad band (0.3-100 keV) spectrum. We found the 
best match with composite models that consist of a combination 
of a number of optically thin plasma (mekal) components plus a 
Gaussian at 6.4 keV, including a total (Wabs) and a partial cov- 
ering absorber (Pcfabs). The statistical significance of the inclu- 
sion of additional mekal components was always verified with a 
F-test. All best-fit spectral parameters are reported in table 4. 
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At least two characteristic temperatures were always 
present, three in the case of IGR J0839, IGRJ1830, IGRJ1650, 
IGRJ1817 and IGRJ1719. We broadly divided these optically 
thin components into a cold (kT c ~ 0.1 - 0.2 keV, exception is 
IGRJ0839 with kT c ~ 0.7 keV), a medium (kT m -4-10 keV) 
and an hot (kTh -30-60 keV) one. The presence of low tem- 
perature plasma is supported by the RGS spectral analysis of 
IGR J1719, IGR J1509, and XSS J0056 (see sect. 3.2.3). 

However, in V2069 Cyg the cold mekal temperature was un- 
constrained, while for RX J0636, the abundance was unreliably 
low (A z -0) A more physically-justified model is obtained by 
substituting the cold mekal with a blackbody (BB) component 
at a temperature of kT B B ~ 70-80 eV (see also figure 9). 

All metal abundances were found to be consistent with the 
solar value except IGR J1 8 17, for which a highly sub-solar 
(0.20 ± 0.07) value was found. 

For all sources we found significant presence of an emission 
feature at about 6.4 keV (see figure 9), with no significant shift 
in energy. Hence, we fixed the centroid energy of the Gaussian 
component at 6.4 keV. The equivalent widths (EW) of this emis- 
sion feature were found -1 00-220 eV, consistent with fluores- 
cent Fe K„ line. This suggests the presence of a continuum re- 
flected component, whose inclusion (reflect), however, pro- 
vides meaningful fit only for IGR J1 650, lowering the tempera- 
ture of the hot mekal (see last line of table 4). 

All sources are affected by strong complex absorption. The 
hydrogen column density (Nh) of the total absorber was found 
to be lower or consistent with the interstellar value in the di- 
rection of the sources. The partial absorber, with typical cov- 
ering fraction - 60%, was found at much higher column den- 
sities (up to 10 23 cm -2 ) and hence is instead intrinsic to the 
sources. For IGRJ1817 and XSS J0056, which were found with 
the highest value of intrinsic absorption (50 ± 10 x 10 22 and 
300 ± 200 x 10 22 cm -2 ) we used the PwabXSPEC model 
(Done & Magdziarz 1998), that assumes a power-law distribu- 
tion of neutral absorbers and is likely more appropriate in case 
of high density intervening absorbing columns. This provided a 
better fit in terms oix 1 - However, the column density resulted to 
be poorly constrained in XSS J0056. 

IGR J0839 is the only source which, in addition to the two in- 
tervening absorbing columns, also displays an evident absorp- 
tion feature found at 0.76 keV. We modeled it using an edge 
component, which statistical significance, evaluated with an 
F-Test, resulted to be 5.6cr. This feature is consistent with being 
an OVII K-shell absorption edge with an optical depth r=0.6. It 
indicates an ionized absorber along the line of sight. 


IGR J 1 7 1 9—4 1 00 
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Fig. 9. Broad-band unfolded spectra, with post fit residuals in the 
lower panels, for three representative cases. Spectra have been 
rebinned for plotting purpose. Composite model (solid black 
line), individual components (dotted blue lines). Up: IGR J1719. 
Wabs*Pcf a bs*(3mekal+gauss) model used. Central: V2069 Cyg. 
Wabs*Pcf a bs*(bbodyrad+mekal+gauss) model used. Bottom: 
IGR J0839. Wabs*Pcfabs*(3mekal+gauss) model used. An ab- 
sorption feature at -0.76 keV is present in the residuals (lower 
panel). 
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Table 4. Spectral properties for the best fitting model. 
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Table 5. column 2\ WD masses as derived from maximum post- 
shock temperature, using Suleimanov et al. (2005) model (see 
text), column 3: Lower limits to the mass accretion rate adopting 
the source distances reported in column 4 and the orbital periods 
in the last column of table 2. 


Source 

M wd 

M q 

M 

Mo yr -1 

D* 

pc 

IGR J0839 

0.95+0.08 

3. 2x10"'° 

1200 

IGR J 1830 

0.85+0.06 

2.4x10-" 

100 

IGR J 1650 

0.92+0.06 

4.1x10-" 

270 

IGRJ1817 

0.96+0.05 

2.4x10-'° 

610 

IGRJ1719 

0.86+0.06 

1.9x10-" 

130 

V2069 Cyg 

0.82+0.08 

5.3x10-'° 

930 

RX J0636 

0.74+0.06 

5.4x10-'° 

1100 

IRG J1509 

0.89+0.08 

2.5x10-'° 

600 

XSS J0056 

0.79+0.07 

2.5x10-'° 

320 


* Lower limits to the distances are determined as follows: we used the 
lower uncertainty values of the adopted orbital periods and the 2MASS 
K band (reddening correction not applied) magnitudes of the optical 
counterparts. For the secondary stars we used the K band absolute mag- 
nitudes predicted for CV donors as a function of P or t, (Knigge 2006), 
assuming that these donors totally contribute to the K band flux. 


Since the hot temperature should be regarded as a lower limit 
to the maximum temperature of the PSR, we also used the model 
of Suleimanov et al. (2005), that takes into account the growth of 
pressure toward the WD surface and hence the change of gravity, 
to obtain a more reliable value for the maximum temperature and 
consequently an estimate of the WD mass. This model is com- 
puted for the continuum only. We included in the fit the same 
absorber components as in the multi-temperature fit, as well as 
a broad Gaussian that takes into account the iron line complex. 
We obtain in all cases an acceptable ^ ed between 1.1-1. 3. The 
resulting masses are found in the range 0.74-0.96 M e (see table 
5). We also note that the temperature of the hot component (and 
hence the masses) determined with the multi-temperature model 
are within 3cr consistent with those found with the PSR model, 
except for XSS J0056 whose temperature would give a unreal- 
istic low WD mass. This could suggest that the hot component 
is close to that at the shock. All other best fit parameters are 
compatible, within uncertainty, with that reported in table 4. For 
IGR J 1 650, we also tried to include a reflection component with- 
out any fit improvement. This is likely due to the lack of cover- 
age in the 10-20keV range and the low statistics above 20keV, 
that prevent us from drawing conclusions on the reflected con- 
tinuum. 

3.2.2. Pulse Phase Spectroscopy 

To understand the role of spectral parameters in producing the 
observed X-ray pulses, we carried out a phase resolved spectral 
analysis. For all sources we extracted the spectra at the maxi- 
mum and minimum of the pulsation (see table 6 for details on 
the interval selection). We kept fixed the column density of the 
total absorber, the metal abundance, and the parameters of the 
hot mekal component to the values found from the phase av- 
eraged spectral fits. We also fixed the intermediate mekal tem- 
perature, since it was always consistent with that derived from 
average fit (except for IGRJ1650). Furthermore, inspection of 
the fluorescent Fe line for phase-dependent velocity shifts gave 
no significant change and hence we fixed the Gaussian centroid 


at 6.4 keV. All other parameters were left free to vary (see table 

4). 

In all cases the spectrum at the minimum is more absorbed 
than that at the maximum (see table 6), due to higher values of 
the hydrogen column density and of the covering fraction (how- 
ever for XSS J0056 we fixed the density of the absorber to the 
average spectrum values). No variation in the normalization of 
the cold mekal component is found except for IGR J1 650 where 
it decreases at minimum. We then conclude that the pulses are 
mainly due to phase-dependent absorption. 

Interesting is the case of V2069 Cyg that instead shows an 
increase in the normalization of the BB component at pulse min- 
imum, in accordance with the behaviour of the energy dependent 
pulses (figure 6). This suggests that the secondary maximum at 
0-0.8, best seen below 1 keV, is produced by an increase of the 
normalization of this optically thick component. 

We found evidence for phase variation in the EWs of the 
Fe fluorescent line only in IGRJ0839 and IGR J1 509. In these 
sources the EWs change between 0.08 ±0.03 and 0. 1 6±0.03 keV 
and 0.09 ± 0.03 and 0.23 ± 0.04 keV at maximum and minimum, 
respectively. 

3.2.3. RGS spectrum 

The RGS spectra of IGRJ1719, IGRJ1509 and XSS J0056 (see 
figure 9) were fitted in two ways. In one, we used continuum 
only models (absorbed BB plus bremsstrahlung) and fitted the 
strongest lines using Gaussians, to derive their centroid, width, 
and intensity in a manner that does not strongly depend on 
the choice of models. In IGR J1 7 19, the lines of OVII, OVIII, 
NeIX and NeX were all strong enough for such a treatment; in 
IGR J1 509, only the two oxygen lines; and for XSS J0056 only 
the OVII line (see table 7). In addition to the line fluxes, which 
are consistent with the prediction of the EPIC model, we found, 
for IGRJ1719 and IGR J1 509, that the lines are mildly broad- 
ened with Gaussian cr of order 1000 km s" 1 (or FWHM ~2500 
km s -1 ) and that the centroids of the He-like triplets (OVII and 
NeIX for the fonner and OVII only for the latter) are between the 
expected energies of the resonant and intercombination compo- 
nents, when fitted with a single Gaussian. 

We also explored global fits, based on the best-fit model for 
the EPIC spectra. For XSS J0056, there is a marginal indica- 
tion that an addition of a OVII edge would improve the fit, even 
though this was not statistically significant in the EPIC fit. This 
should be explored further using future, higher signal-to-noise 
spectra. For IGR J1 7 19, although the basic features of the RGS 
spectra were reproduced using the EPIC model (including the 
normalization), the reduced^ 2 is 1.32 for 363 degrees of free- 
dom. We obtained an improved fit (x 2 = 1.03, 351 do!) by in- 
cluding a Gaussian deconvolution using the XSPEC gsmooth 
model and by letting the mekal density parameter be fitted. The 
fit converged to a high density regime (i.e., with no forbidden 
components for the He-like triplets) with cr ~ 1250 km s -1 , 
consistent with our findings from the phenomenological model. 
Similarly, a high density, mildly broadened model was more suc- 
cessful in fitting the RGS spectra of IGR J1 509. 


4. Discussion 

The CVs analyzed here are all found to be strong X-ray pulsators 
with periods ranging from -470 s to -1800 s. In IGRJ0839, 
IGR J1 8 17, IGR J1 7 19, V2069Cyg, IGR J1 509 and XSS J0056, 
periodicities are also identified in simultaneous OM optical (B 
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Table 6. Spectral parameters of max and min spectra. Other parameters are fixed to their best-fit values of the average spectrum. 


source 

N ffpc 

10 22 cm" 2 

cvf 

% 

kT c 

keV 

n c 

10“ 3 

Am 

10" 3 

n h 

10“ 3 

EW 

keV 

Fo.3-15 

lO" 11 

T 2 /dof 

IGR J0839" 

Max (</> = 0.0 - 0.4) 

6(1) 

40(3) 

0.69(fix) 

0.34(7) 

1.2(3) 

3.9(3) 

0.08(3) 

0.92(3) 

1.20/485 

Min (</ = 0.5 - 0.9) 

7(1) 

81(1) 

0.69(fix) 

< 0.4 

0.7(3) 

3.1(2) 

0.16(3) 

0.60(3) 


IGR J 1830 

Max (</> = 0.9- 1.1) 

5(2) 

49(3) 

0.21(1) 

5(1) 

1.9(3) 

4.7(3) 

0.22(4) 

1.11(6) 

1.03/131 

Min (</ = 0.35 - 0.65) 

14(2) 

65(2) 

0.16(1) 

30(10) 

1.3(2) 

3.9(5) 

0.23(8) 

0.76(4) 


IGR J 1650 

Max (<p = 0.8- 1.1) 

13(2) 

56(2) 

0.09(1) ; 0.7(1)* 

2.8(4) 

0.10(5) 

8.6(5) 

0.25(5) 

1.58(4) 

0.93/353 

Min (<fi = 0.2 - 0.5) 

12(1) 

66(1) 

0.15(4) ; 8(2)* 

0.6(2) 

1.8(5) 

6.6(7) 

0.22(3) 

1.44(6) 


IGR J 1 8 1 7 

Max {</> = 0.55 - 0.75) 

50(15)** 

-0.25(8)^ 

0.26(fix) 

9.6(3) 

15.2(6) 

10(4) 

0.15(4) 

1.7(1) 

1.06/356 

Min (</ = 0.0 - 0.2) 

80(20)** 

-o.i(i/ 

0.26(fix) 

16(7) 

13(9) 

9(3) 

0.10(3) 

1.4(1) 


IGRJ1719 

Max (c/ = 0.2 - 0.5) 

6.7(7) 

39(2) 

0.18(1) 

1.2(1) 

6(1) 

19(1) 

0.14(2) 

4.75(3) 

1.03/1768 

Min (<f> = 0.7- 1.0) 

7.9(6) 

47(1) 

0.15(1) 

1.7(4) 

5(1) 

20(1) 

0.13(2) 

4.64(3) 


V2069 Cyg 

Max(</> = 0.1 -0.55) 

5.4(6) 

61(7) 

0.073(1)** 

5.8(6/ 

- 

7.2(1) 

0.24(8) 

1.33(3) 

1.05/275 

Min (<p = 0.6 - 0.95) 

11(1) 

70(1) 

0.065(1)** 

13(1/ 

- 

6.6(1) 

0.27(4) 

1.07(4) 


RX J0636 

Max (p = 0.3 - 0.7) 

9(2) 

49(3) 

0.08(3)** 

0.6(1/ 

- 

5.2(4) 

0.08(3) 

1.04(3) 

1.08/441 

Min (</> = 0.8 - 1.1) 

8(2) 

64(3) 

0.08(1)** 

0.7(3/ 

- 

4.7(3) 

0.07(6) 

0.87(5) 


IGR J 1509 

Max (0 = 0.75- 1.1) 

5(D 

27(2) 

0.11(1) 

3.2(6) 

- 

11.6(3) 

0.09(3) 

2.50(4) 

1.11/824 

Min (</ = 0.2 - 0.6) 

9(1) 

67(1) 

0.12(1) 

2.8(6) 

- 

10.7(3) 

0.23(4) 

1.82(2) 


XSS J0056 

Max (<f> = 0.4 - 0.7) 

300(fix)* 

-0.64(1)*' 

0.16(2) 

8(2) 

50.0(2) 

- 

0.12(3) 

3.2(1) 

0.96/360 

Min (</ = 0.0 - 0.2) 

300(fix)* 

-0.52(1)*' 

0.12(4) 

20(10) 

48.6(2) 

- 

0.15(3) 

26(1) 



** Max: E erfge = 0.75 ± 0.01 r = 1.4 ± 0.2; Min: E edge = 0.76 ± 0.01 r = 0.9 + 0.2. 
h Medium mekal temperature. 

d P- 

** Blackbodyrad temperature. 

f Blackbodyrad Normalization xlO 4 . Blackbodyrad normalization is defined as Rj; m /D 2 0 , where Dio is the source distance in unit of 10 kpc. 


or U bands) or UV photometry. X-ray pulses are unambiguous 
signature of magnetically confined accretion and hence of the 
WD rotational period. Therefore, all sources can be classified as 
CVs of the IP type. 

The XMM-Newton data have allowed us: to identify for 
the first time the WD spin period in IGRJ1830, IGRJ1650, 
IGR J 1 8 1 7, IGRJ17 19, V2069 Cyg and RX J063 6, and to refine 
the spin period of IGR J0839 and IGR J1 509 (but we determined 
a slightly different period for XSSJ0056). Furthermore, we 
could infer the orbital periods (except for IGRJ1650). We pro- 
vide first estimates for IGRJ0839, IGR J1 830, IGRJ1817 and 
find consistency with the spectroscopic periods for IGR J1 7 19, 
IGR J1509, V2069Cyg and XSS J0056 except RXJ0636. 

We have also presented the first broad band X-ray spectral 
analysis of these sources as well as spin phase-resolved spectra, 
allowing us to identify commonalities to IPs, but also peculiari- 
ties in some sources. 

We discuss here the results in terms of accretion and emis- 
sion properties. 

4. 1. X-ray pulses as tracers of accretion mode 

The presence in the power spectra of a strong and dominant sig- 
nal at the WD spin period, especially in the hard (E>4 keV) X- 


rays, is a signature that accretion occurs predominantly through a 
disc because the circulating material loses memory of the orbital 
motion (Wynn & King 1992). For IGRJ0839, V2069Cyg and 
IGR J 1 6500, we have detected only the spin frequency in X-rays, 
and then they likely are pure disc-accretors. Noteworthy, the OM 
optical (but not the UV) data of IGR J0839 and the ground-based 
optical photometry of IGR J1 6500, IGR J1 71 9 (Pretorius 2009) 
and RX J0636 (Gansicke et al. 2005) reveal instead stronger sig- 
nals at the orbital sideband, cj-Q, rather than at the spin period. 
Since the optical light originates from reprocessing of X-rays at 
fixed regions within the binary (Warner 1995), this makes the 
X-rays the best tracers of the true WD rotational period. 

Six out of nine systems, while displaying a dominant spin 
periodicity in the X-ray band, also show signals at sideband fre- 
quencies (table 2), with spin-to-sideband amplitude ratios ~ 1 .4- 
2.6. Hence, although the bulk of accretion onto the WD takes 
place through a disc, a fraction of it (~35% up to ~70%) passes 
over the disc and is directly channeled onto the WD magnetic 
poles. This hybrid mode of accretion, the disc-overflow, is also 
observed in other confirmed IPs (Hellier 1995; Norton & Taylor 
1996; Norton et al. 1997) and in more recently identified mem- 
bers (Anzolin et al. 2008). Worth mentioning is the case of 
IGR J1 71 9 whose X-ray light curve reveals spin and sideband 
pulses only during the first 2/3 of observation, while no strong 
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Table 7. Spectral parameters with RGS. 


Line 

E 

IGRJ1719 

cr 

Flux 


keV 

eV 

photons cm -2 s _l 


OVII 

0 5706+ 0001 
' j,J/uu - 0 . 0008 

3.8 + 0.7 

(3-3±o 6 ) x 10 -3 

OVI1I 

NeIX 

0 6516+ 00018 

0-9209+“ 

4 3+ 1 - 6 
1*4 

11 3+ 8 - 5 

(7.7±j|) x 10 -4 
(3.7±+V 8 ) x 10“ 4 

NeX 

1 0228±“S 

< 32.2 

(1.3±+i°)xl0- 4 

= 

1250±^ [km/s] 




IGR J 1509 


Line 

E 

cr 

Flux 


keV 

eV 

photons cm -2 s -1 

OVII 

OVIII 

0 5705+ 0 0021 

/WJ -Q.QQ15 

0 6549+™ ^ 

4-l ± u 

2-8±iii 

(7.3±2 ®) X lO” 5 
(3.6±‘|)x 10- 5 

°g = 

1240±5qq [km/s] 




XSS0056 


Line 

E 

cr 

Flux 


keV 

eV 

photons cm -2 s _l 

OVII 

0 5676+ 00618 
U,JU/UX 0 . 0444 

< 240 

(3.7±J° 2 4 )x 10“ 4 

coherent signals are found afterward. This system is the one 

a low spin to beat amplitude ratio of ~ 

1.4 and a spin PF~ 1 


It could be possible that the two variabilities interfere destruc- 
tively, as also suggested by the large jumps in the phase of the 
signal with time (see figure 2) and the low-frequency variations 
in the power spectrum. Hence, the somewhat peculiar behaviour 
of this IP could be due to substantial disc-overflow. 

Two systems, IGR J1817 and XSS J0056, reveal a weak sig- 
nal at a putative positive sideband, at + Q, that would imply a 
retrograde WD rotation. We cannot exclude that the true spin 
period is the latter and hence the stronger signal is the beat. In 
this case the two systems would be beat-dominated at both X- 
ray and optical wavelengths and then disc-less accretor. A simi- 
lar problem raised for BG CMi (Norton et al. 1 992; Garlick et al. 
1994; de Martino et al. 1995). To solve this issue, both fast time- 
resolved spectroscopy and polarimetry are needed. 

Furthermore, for IGRJ1830 and RXJ0636 the presence of 
the a) - 20. sideband would imply an orbital modulation of the 
amplitude of the beat frequency if this is an intrinsic emission. 
Hence, the beat and the spin pulses would be produced by dif- 
ferent modes, stream-fed and disc-fed accretion respectively, see 
also Norton et al. (1997). 

Long-term variation consistent with the orbital period was 
found in XSS J0056, and likely in V2069Cyg. In the latter, we 
also detected a dip at minimum of the possible orbital modu- 
lation, more pronounced in the soft bands, indicating a spectral 
hardening. Orbital variations are often found in IPs and are likely 
due to absorption effects from fixed regions within the binary 
frame, such as the disc edge, in systems viewed at relatively large 
inclination angles, i<; 60° (Parker et al. 2005). 

Using the adopted orbital periods, reported in table 2, the 
spin-orbit period ratios, defined as P w /Pn, are lower than 0.1, 



Fig. 10. From top to bottom, IGRJ1719, IGRJ1509, and 
XSS J0056. RGS1 data and model (black), RGS2 data and model 
(blue). The models use the EPIC best-fit values. Those for 
IGR J1 71 9 and IGR J 1509 also include fitting for gsmooth cr 
and for mekal density parameters. Top panel inset shows the 
OVII line wavelength range (20.6-23.1 A) with the EPIC best- 
fit low density mekal (black), the high density, narrow mekal 
(blue) and the high density, gsmooth broadened meakl (red). For 
XSS J0056 we show the EPIC best-fit model with (black) and 
without (red) the edge. 


except for IGR J1830 (0.12), and cluster around 0.03-0.05 (see 
figure 1 1). All but one (XSS J0056), have Pci above the orbital 
gap. A strong asynchronism (P^/Pn < 0.1) and long orbital pe- 
riods are predicted by evolutionary models, where mCVs evolve 
towards shorter orbital periods and, depending on the WD mag- 
netic field strength, may or may not eventually synchronize 
(Norton et al. 2004, 2008). In long period highly asynchronous 
systems, the accretion flow is expected to take part in a disc-like 
configuration, as indeed inferred from X-rays. Such configura- 
tion may change while the system evolves towards short orbital 
periods, depending on the WD magnetic field strength. To fur- 
ther progress in this issue, polarimetric measures of new systems 
are crucial to assess whether these systems will eventually syn- 
chronize. For IGR J 1 509 a recent polarimetric survey of new IPs 
has indeed revealed that the optical-red light is polarized at a few 
percent (~ 2%), suggesting a strong magnetic field (>10MG) 
and a likely evolution towards synchronism (Potter et al. 2012). 

4.2. Origin of pulses 

In most of the sampled sources the harmonics of the fundamental 
(up to 4a>) are detected, indicating departures from strictly sinu- 
soidal pulse shapes. Non strictly sinusoidal pulses are seen also 
in the hard X-rays, where absorption is negligible. This suggests 
that the emitting regions depart from a symmetric dipole struc- 
ture. 
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V2069 Cyg is the one showing the most structured pulse in 
our sample (see figure 6), with a secondary maximum centered 
on a primary broad minimum seen at all energies. This indi- 
cates the presence of a secondary emitting pole. Furthermore, 
FIRs in the 0.3-1 keV and 1-3 keV bands display a softening 
at spin minimum, indicating an additional contribution at this 
phase. V2069 Cyg also shows a BB component at soft energies 
(see Sect. 3.2.1) whose contribution increases at the minimum 
of pulse (see table 6). Hence, a reprocessed component at the 
WD surface is present and is anti-phased with the optically thin 
emission from the PSR. 

On the other hand, IGRJ1817 pulsates at a period that is 
half that found at optical wavelengths, indicating two symmet- 
ric poles contributing almost equally in the X-rays. The optical 
pulse is broad and overlaps the two X-ray maxima (figure 6) and 
hence likely originates from wide upper accretion regions above 
the two poles. 

For all sources the pulses are energy dependent, with ampli- 
tude decreasing at high energies. The softest band (0.3-1 keV) 
has pulse fractions even up to 99%, seen in IGR J0839, decreas- 
ing to a few percent in the harder band 5-15keV. The pulses 
harden at spin minimum in most sources as confirmed from the 
spectral analysis (see Sect. 3.2.2). This is a typical behaviour ob- 
served in IP type CVs and consistent with the accretion curtain 
scenario proposed by Rosen et al. (1988), where the modulation 
is mainly caused by spin-dependent photoelectric absorption in 
an arc-shaped accretion curtain. The absorption is expected to 
be larger when viewing the curtain along the magnetic field lines 
(spin minimum). At these phases the softest regions are absorbed 
resulting in a spectral hardening. 

The curtain is also expected to contribute in the UV and op- 
tical light and hence the pulses at these wavelengths should be 
in phase with the X-rays. This is the case of IGRJ0839, and 
IGR J 1509, though the optical light is dominated by the beat in 
the first system. A peculiar case is XSS J0056 that reveals opti- 
cal pulses changing phase with time: they are phase aligned with 
the X-rays only during the first 2/5th of observation. This may 
suggest a non-stationary contribution from two poles in the op- 
tical light. Also IGR J1 71 9 and V2069 Cyg display anti-phased 
X-ray and optical pulses that could be due either to a more lumi- 
nous secondary pole or to the heated WD polar cap. For V2069 
Cyg the second option could be favored because of the detection 
of the BB component in its spectrum. 

4.3. X-ray emission components and energy budget 

The analysis of the average 0.3-100 keV spectra gives evi- 
dence of multi-temperature emission components as well as of 
complex intrinsic absorption. A low temperature (~0.1-0.7keV) 
plasma traces the regions closer to the WD surface. The oxygen 
H- and He-like lines detected in the RGS spectra of IGR J1 7 19 
and IGR J1 509 give further evidence for the presence of low- 
temperature plasma. Their intensity ratio (Mauche 2002) in- 
dicates a temperature of ~0.2keV and ~0.1 keV, respectively. 
Since no forbidden component is detected in the He-like triplet, 
these lines form in high density regime (n > 10 13 cm -3 ). The 
presence of Ne IX and Ne X and their ratio further indicate a 
temperature at ~0.3 keV in IGR J1 7 19. 

The intermediate temperature component is required in the 
PN spectral fits to account for the K„ lines of Fe XXV (6.7 keV) 
and Fe XXVI (6.97 keV). These are strong in the spectrum of 
XSS J0056 thus providing a satisfactory fit with only two mekal 
components. 


On the other hand, the high temperature component should 
originate in regions close to the shock. The wide range of values 
(~30-60keV) found in our sample would also imply a range of 
WD masses. The emitting hot regions have typically larger (~3 
times) emission measures (EM, normalizations in table 4) than 
those of the intermediate component, being that of IGR J1 650 ~5 
times. As EM oc N 3 l 3 , the larger EMi lot might indicate a larger 
emitting volume and hence, on average, a larger size. 

Hence, the multi-temperature spectra indicate a temperature 
gradient in the PSR though we are unable to map it with the 
present data. 

Different are the cases of RX J0636 and V2069 Cyg that dis- 
play a non-negligible soft BB component with kT B B ~ 70 - 80 
eV. These two sources add themselves to the small group of soft 
X-ray IPs (Anzolin et al. 2008) increasing the roster to 15 sys- 
tems. The ratio of bolometric fluxes L B B/Lth. are ~ 0.85 for 
V2069 Cyg and ~ 0.13 for RXJ0636, respectively. Therefore, 
this component represents a substantial fraction of the energy 
budget in V2069 Cyg, but not in RX J0636. The normalization 
of this optically thick component increases at spin minimum in 
both sources. It can be identified with the X-ray heated WD po- 
lar cap, whose projected area is larger at this phase. We estimate 
an average spot area in V2069 Cyg, a BB = 4.5 x 10 14 D 3 00pc and 
hence a fractional WD area / ~ 7.6 x 10 -5 Dg 00pc , for M WD = 
0.82 (see table 5 and also its notes for the lower limit estimate of 
distances). For RXJ0636 we estimate a BB = 6.3xl0 13 D 3 kpc and 

/ ~ 9 x 10 ~ 6 D 2 lkpc . These values are small, but comparable to 
what found in other soft X-ray IPs with similar BB temperatures. 

An additional component, ubiquitous in magnetic CVs, is the 
K„ fluorescent line at 6.4 keV. It is detected in all sources of our 
sample with large EWs (~ 100-200 eV), arising from reflection of 
cold material. The lack of velocity shifts and the EW variability 
at pulse maximum and minimum in IGRJ0839 and IGR J1509 
suggest a likely origin at the WD surface. Such feature should 
be accompanied by a reflection continuum but our data do not 
allow to constrain this component. 

A lower limit to the mass accretion rate can be derived as- 
suming that the accretion luminosity is totally emitted in the X- 
rays: L accr = GM 0 M WD R^ D > L BB + L th , where L th is the 
bolometric luminosity of the X-ray optically thin components 
and L bb is the reprocessed X-ray luminosity, that is only in- 
cluded for V2069 Cyg and RX J0636. The estimated values are 
reported in column 3 of table 5. The secular mass accretion rates 
expected from angular momentum loss (AML) due to magnetic 
braking at the adopted orbital periods, would be in the range 
of KR V - 10 -9 M 0 (McDermott & Taam 1989), for the sources 
above the gap. Hence, unless they are located at a distance 10 
times larger than the estimated lower limits (very unlikely pos- 
sibly except for IGR J1 830 and IGR J1 7 19), the X-ray emission 
can not be solely representative of accretion luminosity. A sub- 
stantial fraction should be reprocessed and emitted in low en- 
ergy bands, as supported by the detection of periodic signals at 
the spin and/or beat periods at UV and/or optical wavelengths, 
as firstly raised by Mukaietal. (1994). The only exception is 
XSS J0056 that falls in the period gap, for which the lower limit 
to the rate is compatible with the values expected either from 
magnetic braking or gravitational radiation (Warner 1995). 

The spectral analysis has also shown the presence of dense 
(column densities up to ~ 10 23 cm -2 ) material partially ab- 
sorbing the X-ray emission. The spin phase dependence of ab- 
sorption indicates that the spin pulsations are mainly due to 
photoelectric absorption from pre-shock material in the accre- 
tion curtain. Furthermore, we found indication for the pres- 
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ence of absorbing ionized material in one source (IGRJ0839), 
whose signature is an OVII absorption edge detected at 0.76 keV. 
IGR J0839 adds as the third IP showing this feature (Mukai et al. 
2001; de Martino et al. 2008). A warm absorber is commonly 
recognized in Low Mass X-ray Binaries (LMXB) through the 
presence of absorption edges and lines that have been related to 
ionized atmosphere above the accretion disc, giving rise to dips 
in close to edge-on systems (Diaz Trigo et al. 2006). The spectra 
at pulse maximum and minimum of this source do not provide 
strong evidence of changing in the edge parameters, suggesting 
that this material could be located within the binary, but not in 
the pre-shock magnetically confined accretion flow. A possible 
similarity to the LMXB dippers should be further investigated 
with phase-resolved spectroscopy along the binary period that 
we are unable to perform with the present data. 

4.4. Are hard X-ray IPs different? 

The hard X-ray sources analyzed here are confirmed as IP type 
CVs. Their spin periods are well within the range of previously 
known members of this group as depicted in figure 1 1 where 
our sample is shown together with the current confirmed IPs. 
Orbital periods have been estimated combining X-ray and opti- 
cal results. Most of them are found above the 2-3 hr period gap 
except for XSS J0056 whose high asynchronism (P w /Pn=0.05) 
might suggest it will never synchronize. The inferred spin-to- 
orbit period ratios are £5 0.12 consistent with most IPs. Hence, 
there is no striking difference between hard X-ray discovered IPs 
and other known IPs. Nonetheless, the orbital period distribution 
is now being also populated by the new hard systems at Pn S 5 
h. 

To inspect whether these selected INTEGRAL and Swift IPs 
possess massive WD we derived their masses using the maxi- 
mum temperature obtained from fitting the broad-band spectra 
with the PSR model and taking it as an estimate of the shock 
temperature (see table 5). These are obtained using T s h 0C k = 
3/8GM WD /u m H/kR WD and adopting the WD mass-radius rela- 
tion from Nauenberg ( 1 972). The mean WD mass of our sample 
is 0.86±0.07 M 0 . Consistent results were also found in a recent 
study of 17 IPs made by Yuasa et al. (2010), though they sug- 
gest a likely bias toward high masses (due to a high energy sam- 
pling). The average WD mass we found is also fully consistent 
with a recent determination of CV primary masses of 0.8±0.2 
made by Zorotovic et al. (201 1). These authors also find that the 
WDs in CVs are substantially more massive than those in pre- 
CVs (0.67±0.21M o ) and than single WDs (^0.6 M 0 ). While the 
latter cases imply that the WDs increase their mass during the 
CV stage, the fact that IPs have primaries with masses similar to 
those of the whole CV population, suggests that the mass is not 
the driving parameter. Therefore, we conclude that the hard X- 
ray detected CVs, though possibly biased toward high mass, do 
not show remarkable difference from the whole CV population. 

With the increase in the INTEGRAL and Swift exposures 
with time, most of the undiscovered IPs within ~1.5 kpc will 
be likely detected. This is the case of the faint RXJ0636, 
which only recently appeared in the new Swift/BAT catalogue 
(Cusumano et al. 2010). Therefore, the true population of IP 
type CVs is expected to be eventually unveiled by sensitive hard 
X-ray surveys of next generation mission. Moreover, as the new 
discovered IPs are being found in the previously poorly popu- 
lated range of long orbital periods (P>5 h), this will allow us 
a better understanding of the evolution these systems through 
comparison with recent models, e.g.Knigge et al. (2011). 


5. Conclusions 

The main goal of this work is the characterization of a sample 
of 9 new hard X-ray selected CVs, to unambiguously identify 
them as magnetic systems of the IP type. The main results can 
be summarized as follow: 

- All sources are strong X-ray pulsators with periods between 
470-1 820 s and thus are confirmed CVs of the IP type. 

- All but two sources are spin dominated systems in the X- 
rays. IGR J1 8 17 and XSSJ0056, remain ambiguous cases. 
IGRJ0839, IGRJ1719, IGRJ1650 and RXJ0636 are beat 
dominated at optical wavelengths. 

- IGRJ0839, V2069Cyg and IGR J1 6500 are disc-fed accre- 
tors, while IGR J1830, IGR J1719, RX J0636 and IGR J1509 
display a disc-overflow accretion mode. IGR J1 8 17 and 
XSS J0056 could be either stream-fed or disc-fed systems. 

- We also estimated the orbital periods (except for IGR J1 650), 
all found to be above the 2-3 h period gap. XSS J0056 is 
the only exception, lying in the gap. We derived first esti- 
mates for IGR J0839 (8 ± 1 h), IGR J1 830 (4.2 ± 0.2 h), and 
IGRJ1817 (6. 3-8. 7 h). 

- The amplitude of the X-ray pulse was always found to de- 
crease with energy due to complex absorption. 

- We identified two or three optically thin components with 
characteristic temperatures in the ranges: 0. 1-0.7, 4-10 and 
30-60 keV (signature of a temperature gradient in the post- 
shock region). Three are present in six sources. 

- V2069 Cyg and RX J0636 are instead found to possess a soft 
X-ray optically thick component at kT BB ~80 eV. This in- 
creases to 15 the current roster of soft X-ray IPs, confirming 
a relatively large (~ 30%) incidence. 

- A strong K ff Fe line at 6.4 keV is present in all sources likely 
originating at the WD surface from the reflection of X-rays 
from cold material. 

- The spectrum of IGRJ0839 reveals an absorption edge at 
0.76 keV from OVII. This increases to three the IPs were a 
warm absorber is found. 

- The WD masses of our sample are found between 0.74 and 
0.96 M 0 , with an average mass 0.86 ± 0.07 M 0 , suggesting 
no striking difference with the whole CV population. 
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